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Abstract We review the anomalous behavior of solid He at low temperatures with partic- 
ular attention to the role of structural defects present in solid. The discussion centers around 
the possible role of two level systems and structural glassy components for inducing the 
observed anomalies. We propose that the origin of glassy behavior is due to the dynamics of 
defects like dislocations formed in 4 He. Within the developed framework of glassy compo- 
nents in a solid, we give a summary of the results and predictions for the effects that cover 
the mechanical, thermodynamic, viscoelastic, and electro-elastic contributions of the glassy 
response of solid 4 He. Our proposed glass model for solid 4 He has several implications: (1) 
The anomalous properties of 4 He can be accounted for by allowing defects to freeze out 
at lowest temperatures. The dynamics of solid 4 He is governed by glasslike (glassy) relax- 
ation processes and the distribution of relaxation times varies significantly between differ- 
ent torsional oscillator, shear modulus, and dielectric function experiments. (2) Any defect 
freeze-out will be accompanied by thermodynamic signatures consistent with entropy con- 
tributions from defects. It follows that such entropy contribution is much smaller than the 
required superfluid fraction, yet it is sufficient to account for excess entropy at lowest tem- 
peratures. (3) We predict a Cole-Cole type relation between the real and imaginary part of 
the response functions for rotational and planar shear that is occurring due to the dynam- 
ics of defects. Similar results apply for other response functions. (4) Using the framework 
of glassy dynamics, we predict low-frequency yet to be measured electro-elastic features 
in defect rich 4 He crystals. These predictions allow one to directly test the ideas and very 
presence of glassy contributions in 4 He. 
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1 Introduction 

The discovery of anomalous frequency and dissipation behavior seen in torsional oscilla- 
tors (TOs) GlU at low temperatures in 4 He has stimulated numerous investigations. These 
anomalies have been argued to demonstrate a non-classical rotational inertia (NCRI) of 
the l ong ago pre dicted supersolid quantum stated @>& Successive TO experiments 
Il5ll confirmed the finding of the reported anomalous behavior. Hys- 
teresis behavior and long equilibration times have been observed fl(il[l3l[l6ll . which depend 
strongly on growth histor y an d annealing @J. In the same temperature range, experiments 
including shear modulus l!7L fl8ll . ultrasonic lfl9l l20ll and heat propagation fl9ll have also 
shown various anomalies. The character and the existence of mass flow is still a matter of 
intense investigation. Experiments designed to probe for mass flow by squeezing the lattice 
report no such flow Il^[l^l2ll . l23 . l23ll24l . l2^l23l . However, experiments in which a chem- 
ical potential gradient was created via coupling to a superfluid reservoir, suggest mass flow 
of an unusual type l24l I27I I2H , 12^, l30ll. These and many other results have led to a flurry 
of activity. Structural measurements l3ll. l32ll suggest that solid 4 He may be composed of a 
dynamic mosaic of crystals, highlight the importance of defects, and report the absence of 
notable structural change in the vicinity of the putative supersolid transition. The intricate 
structure and dynamics of sold 4 He make it a fascinating system. We note that when a mosaic 
of 4 He crystals with intervening liquid channels is placed in a metallic container 
thermal expansion effects (in particular the larger thermal expansion coefficient of helium 
vis a vis that of the metallic container) may be at play in blocking a remnant superleak. 
Such a thermal compression will effectively shut down any mass flow and close the open 
channels as the temperature is raised. This scenario remains a viable explanation for the 
anomalous mass flow and fountain effect reported by Hallock's group l24l[27ll28l. [2^ . [30ll . 
until a superleak can be ruled out. 

After eight years of intense experimental and theoretical investigations one must ask 
what are the established facts, what are the current expectations and hypotheses, and what 
are the future directions of research in solid 4 He. There are numerous reviews and progress 
updates available describing the field l34l. I35L [3^1 and addressing some of these questions. 
Most of the literature reviewing the subject deals with the notion of supersolidity in 4 He 
as established and proceeds with the discussion on the status and future experiments with 
the goal of further "proving" the existence of a supersolid phase transition at very low tem- 
peratures. Amongst the many exciting proposals concer ning 4 He as well as supersolids, we 
briefly mention Andreev's proposal for superglass l37ll38l [39ll , Anderson's suggestion of 
vortex proliferation and flow 14(1 kill and supersolid dislocation cores I42l.l43l, l44l . l45tl . As 
will become evident in later sections, our analysis centers on the dynamical effects of de- 
fects and as such may include vortices, dislocations, or any other defects. Estimating the 
product of typical dislocation core sizes (~ 3 nm) and their density (~ 10 10 m~ 2 ), it is seen 
that a direct NCRI origin from superfluid dislocation cores in 4 He cannot account for the 
magnitude of the TO anomaly (it is orders of magnitude too small). 

We believe that at this mature stage of research in 4 He a somewhat more general view 
on possible options is needed, where one asks the question: what are the options for possible 
states that might form at low temperatures in 4 He and by implication in other solid bosonic 
matter? By taking a broader view one explicitly allows for states other than pure supersolid- 
ity and includes possible coexistence phases to form as well. This review is a contribution 
in the spirit of broadening the conversation by explicitly allowing for other components or 
"active ingredients" to be present in addition to, or perhaps instead of, supersolidity in 4 He. 
Thus, if one accepts the presence of defects in a solid, then naturally the question arises 
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Fig. 1 The phase diagram of 4 He with the 
putative supersolid phase transition below 
0.3 K d or an alternate order-disorder 
dominated crossover region governed by 
impedance matching between the applied 
frequency and internal relaxation processes. 
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about the dynamic signatures of such crystal defects and whether they can dominate the 
response to an external stimulus. 

We provide a brief critical analysis of some of the existing data and point out that a 
significant fraction of the data on mechanical, thermodynamic and dielectric properties of 
4 He can be analyzed in terms of the emergence of a dissipative viscous component that we 
shall call glassy component l4q. l47l . l48ll49ll5Cll . This extra component by itself is sufficient 
to modify numerous properties of solid 4 He and can be responsible for anomalous thermo- 
dynamic, elastic and dielectric properties in solid 4 He observed in experiments. The glassy 
component is not a supersolid in the classical sense f34ll . yet it can coexist and couple to a 
supersolid component as some of the proposed superglass phases indicate fl3l . l5ll[52 . l53ll . 

The precise nature of the state of solid 4 He at lowest temperatures remains a puzzle. 
In the past, investigation of 4 He has played an important role in the development of basic 
concepts in modern condensed matter physics like superfluidity, order parameter, topolog- 
ical excitations, and critical exponents. Given such a prominent role it played in the past it 
is paramount to come up with the resolution of the puzzles that are clearly seen in experi- 
ments at lowest temperatures. Yet there is one important difference that might be key to a 
solution this time. We propose that precisely because the compound is clean and very well 
characterized the enabling component for the anomalies at lowest temperature are defects 
that undergo freeze-out and constitute a glass like component. 

To sharpen this point of the discussion, any supersolid component would imply some 
sort of two-fluid hydrodynamics that schematically equates the total mass current of helium 
atoms as a sum of normal component and superfluid component determined by the normal 
and superfluid mass fractions p„ and p, and respective velocities. In this two-fluid picture 
the total mass current is given as 

j =P«v„+p s v s (1) 

In any superfluid or supersolid phase the coefficients p n (T,P) and p s (T,P) are functions of 
thermodynamic variables like temperature T and pressure P. There is no direct evidence of 
either dc or ac supersolid mass current at lowest temperatures, where the putative supersolid 
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state sets in I8L \rft\ . except for the experiments by Hallock I24I I27L I28L l3(ill . However, it 
remains to be seen whether Hallock's results of mass flow and fountain effect cannot be 
explained through the presence of superleaks connecting the superfluid leads through solid 
helium with channels. Therefore, some even proposed that one cannot "squeeze a superfluid 
component out of a stone" 1 54] . We know that these supersolid expectations do not apply in 
the case of solid 4 He at lowest temperatures, because TO experiments clearly indicate that 
the period and damping exhibit significant hysteresis effects and strong dependence on the 
history of sample preparations and annealing. All these experimental observations combined 
would suggest to an impartial observer that there is at the very least, in addition to super- 
solidity, another physical component at play in 4 He. Our analysis of various experiments 
suggests that the crossovers, seen in the specific heat, TO, shear modulus and dielectric 
function experiments ITfjl l55ll53 , l57ll are a reflection of the physics whose origin is not due 
to supersolidity, but a consequence of the dynamics of defects in solid 4 He . 

Specifically, we propose the presence of a glassy component in solid 4 He at low tem- 
peratures in order to explain the observed anomalous linear temperature dependence in the 
specific heat of an otherwise perfect Debye solid 4 He J4rj|4^l . We used similar ideas when 
we analyzed the TO, shear modulus, and dielectric properties by assuming the presence of a 
glassy component at the parts-per-million (ppm) concentration and asked what the dynamic 
consequences should be. With the wealth of data available we do not attempt to provide a 
complete overview of the field, but give a summary of the work centered on the role of a 
glassy component in an otherwise nearly perfect crystal. 

The exact nature of the glassy component is not known. For example, it may be caused 
by two-level systems of pinned dislocation lines, vortex excitations, etc. It is however impor- 
tant to point out that the amplitude of period shift can be changed dramatically and depends 
on growth history and annealing procedures of the crystal. To explain the puzzling fea- 
tures of solid 4 He we conjecture that structural defects like localized dislocation segments 
or groups of displaced (out-of-equilibrium) atoms effectively form a set of two level sys- 
tems (TLSs) which are present at low temperatures. These immobile crystalline defects will 
affect the thermodynamics of solid 4 He| 46, the mechanics of TOs 1 47] and shear modu- 
lus J49I . I50I1 . and dielectric properties l58il . Other mechanisms for TLS and glassy dynamics, 
e.g., due to point defects and grain boundaries are possible as well. We are at a stage where 
phenomenology allows us to make progress with testable predictions, while a microscopic 
picture of crystal defects and interactions is still missing. 

Early on it was recognized that pinned vibrating dislocation lines can account for a 
plethora of anelastic damping phenomena in the ultrasound, TO, and shear modulus experi- 
ments. Since most experiments are believed to be in the linear or elastic strain-stress regime 
any plastic deformation due to the motion of gliding dislocation loops has been neglected. 
However, this is not necessarily so. We proposed [59] that dislocation-induced anomalous 
softening of solid 4 He is possible due to the classical motion of gliding dislocation lines in 
slip planes. This picture of dislocation motion is widely accepted in conventional metals. 
Similar effects are at play in the quantum arena where mobile dislocations (dislocation cur- 
rents) lead to a screening of applied shear via a Meissner-Higgs type effect J60l . [6lll . Such 
unpinned dislocations that screen shear render the system more fluid and may, in line with 
the framework that we advance in this article, trigger the 4 He anomalies as the mobile dislo- 
cations become quenched at low temperatures 0,E^]. Recent experiments provide further 
impetus for such a picture I62ll63ll . 

The technique of choice for interrogation of solid 4 He has been the TO with varying 
degrees of complexity. However, the TO does not provide any direct information on the 
microstructure of samples. More direct structural x-ray and neutron measurements do not 
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have the adequate resolution to detect any changes at the ppm level in the structure of 4 He 
at lowest temperatures. Additional challenges arise given how small the volume fraction of 
the glassy component can be. We estimate it to be in the range of few hundreds of ppm 
in the specific heat and pressure contribution. Therefore the precise characterization of the 
microstructure of solid 4 He, growth history, annealing, and 3 He dependence remain pressing 
issues in resolving the hypothesis of the presence of a glassy component and TLS in solid 
4 He. The notion of importance of the role of disorder in solid 4 He received further support 
over the years in observations of the strong dependence of TO results on the specific design 
of sample cells, see work by the group of Chan l64h . It is hard to imagine that an intrinsic 
material property like supersolidity should depend strongly on the stiffness or geometry of 
the torsional oscillator apparatus, while crystalline disorder can easily be affected by those 
design properties. 

In the analysis of the observed excess specific heat, we used a model of independent 
TLS, which gives the canonical signature of a linear in temperature contribution at lowest 
temperatures. Building on the presence of TLS we evaluated the mechanical properties of 
the TO using a model of quenched defects. This approach allows us to make predictions on 
the viscoelastic properties of 4 He and on the electro-elastic coupling that can be tested in a 
setup that does not require the TO and hence can be tested over a much broader frequency 
range. These predictions allow to directly verify the very presence of quenched defects in 
4 He. We therefore would welcome any direct tests of our ideas. 

The remainder of this paper is organized as follows. We present the general discussion 
on the role of quenched defects and disorder in solid 4 He with particular attention to the 
consequence of defects on the dynamics in Section [2] This is followed in Section [3] by an 
analysis of the thermodynamic properties of solid 4 He. In Section |4j we present a unified 
framework to analyze dynamical response function that invokes arbitrarily high order back- 
action effects of defects onto the solid bulk. We then invoke, in Section[5]this approach and 
summarize our analysis of the torsional oscillator. In Section [6] we discuss the shear mod- 
ulus analysis and the strain-stress relations using a viscoelastic model designed to capture 
the anelastic contribution from defects. In Section[7]we discuss predictions for the dielectric 
properties that follow from the viscoelastic model with locally frozen-in defect dynamics. 
Finally we conclude with a discussion and give our view on future perspectives in the field. 

2 Defect Quenching and Its Implications 

In this review, we will examine the general consequences of a transition from mobile defects 
(or dynamic fluid-like components) at high temperatures to quenched immobile defects at 
low temperatures. Microscopically, these defects may be dislocations, grain boundaries, vor- 
tices, or others. It may be posited that in quantum solids such as 4 He , zero point motion 
leads to larger dynamics than is common in classical solids. To put the discussion in per- 
spective, we recall a few rudiments concerning annealing and quenching. When quenched, 
systems fall out of equilibrium. En route to non-equilibrium states, relaxation times increase 
until the dynamic components essentially "freeze" (on pertinent experimental time scales) 
into an amorphous state. In materials with (sufficiently large) external disorder, quenching 
may lead to "spin-glass" characteristics. By contrast, in the absence of imposed external 
disorder, when fluid components (either classical or quantum) fall out of equilibrium by suf- 
ficient rapid cooling (so-called "super-cooling") to low enough temperatures, the resultant 
state is termed a "glass". As liquids are supercooled, their characteristic relaxation times and 
viscosity may increase dramatically. If, instead, the temperature is lowered at a sufficiently 
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slow rate, the system does not quench and instead remains an "annealed" system in ther- 
mal equilibrium. Notwithstanding exciting progress, exactly how the dynamic components 
evolve in 4 He crystals as the temperature is lowered is, currently, an open problem. An ini- 
tially surprising and undeniable feature that has, by now, been seen in many experiments is 
that the solid 4 He anomalies depend dramatically on the growth history of the crystal and 
diminish as the system is cooled down slowly and defect quenching is thwarted. Memory 
effects reminiscent to those in glasses are further present. These observations imply that 
"there is more to life" than static NCRI and other annealed supersolid properties on their 
own- quenching plays a definitive role in triggering the observed effects. To understand the 
experimental observations and build a predictive theoretical framework, we invoke general 
physical principles allowing a computation of response functions and using as input known 
characteristics of quenching. Physically, as alluded to in the Introduction, the quenching 
is that of the mobile defects (which may constitute only a tiny fraction of the system) as 
they become arrested against a crystalline background. Our analysis does not rule out the 
presence of a small supersolid component. The observed large change in the TO dissipation 
cannot be solely described by uniform Bose-Einstein condensation I47l l4^ . l65ll . It remains 
to be seen if nonuniform Bose-Einstein condensation alone either along grain boundaries 
l6r3l or along the axis of screw dislocations I5ll. l67l . l63l can explain the dynamic response 
of TOs. Note that simple estimates of the supersolid fraction of dislocation cores are orders 
of magnitude too small. 

There exists a vast literature on defect quenching in systems that range from vortices in 
superfluid Helium to cosmic strings I69ll70ll and countless others. Our particular focus is, of 
course, on defects in a crystalline system (solid 4 He). Quenched dynamics of such defects is 
associated with a change of plasticity and related internal dissipation. Dielectric (and other) 
response functions in systems of varying plasticity, such as various glass formers as their 
temperature is lowered, indicate, in a nearly universal fashion, the presence of a distribution 
of local relaxation times. These lead to the canonical Cole-Cole or Cole-Davidson distribu- 
tion functions and related forms as we briefly elaborate. In an overdamped dissipative sys- 
tem, an impulse (e.g., an external electric field or an elastic deformation) at time t = leads a 
response g(t) which at later times scales as g s i ng i e ~ exp(— f/f) where z is the (single) relax- 
ation time. When Fourier transformed to the complex frequency (a>) plane, this leads to the 
Debye relaxor gsinglei®) = £o/(l — i(Oz). Now, in systems that exhibit a distribution P(z') 
of local relaxation events, the response functions attain the form / dz'P(z')exp(—t/z'). 
Empirically, in dissipative plastic systems, relaxations scale as (exp(— t/z) c ) with a power 
< c < 1 that leads to a "stretching" (slower decay) of the response function as compared 
to its single overdamped mode form of exp(— t/z). This stretched exponential and other 
similar forms of the response function capture the quintessence of the distribution of re- 
laxation times. Two widely used relaxation time distributions are the Cole-Cole (CC) and 
Davidson-Cole (DC) functions that describe a superposition of overdamped oscillators (De- 
bye relaxors) 17 ll . 17211 . With g((o) = goG(w), where go is a material specific constant, these 
two forms are given by different choices for the function G, 

G cc (co) = \/[\-(icoz) a }, 
G DC {w) = \/[\-mzf. (2) 

Values of a and jS that differ from unity qualitatively play the role of the real-time stretching 
exponent c. In the dc limit the mechanical motion of any mobile component will have ceased 
and there will be no relative motion and no transients. Therefore the coefficient go is gener- 
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Fig. 2 Density of states (DOS) of the two-level tunnel- 
ing system. The black-dashed line represents the DOS 
for the standard glass model (H HI, |4(Jl, while the 
blue-solid line is the truncated DOS used to describing 
the TLS with a cutoff energy. 



standard glass model 
glass model with cutoff 




ally a function of frequency. In the case of the single TO its value is set by the resonance, 
go « go («9rj)- 

These relaxation times can be associated with a distribution of TLSs describing viable 
low temperature configurations of the defects. The simple TLS analysis can account for ther- 
modynamic measurements. Recent work f73h obtained results beyond the TLS model with 
fewer parameters for generic non-uniform systems irrespective of specific microscopic ori- 
gin. We will, for the sake of simplicity, review our work on the low temperature properties of 
4 He assuming TLSs. We conjectured 14^1 that structural defects, e.g., localized dislocation 
segments, form such a set of TLSs at low temperatures. These immobile crystal defects af- 
fect the thermodynamics I46ll48ll of bulk 4 He and the mechanics l47ll of the TO loaded with 
4 He. For the analysis of the specific heat, we used independent TLS to obtain the universal 
signature of a linear-in-temperature specific heat term at low temperatures. 



3 Thermodynamics 

Any true phase transition, including supersolid, is accompanied by a thermodynamic signa- 
ture. Therefore it was anticipated that thermodynamic measurements will resolve the exist- 
ing puzzles of supersolidity. The search for such thermodynami c sig natures proved to be 
challenging so far, see e.g., measurements of the specific heat jl4. fj^L \ldi l55l Wk. r78ll . 
measurements of the pressure dependence of the melting curve I79l . 18011 , and pressure- 
temperature measurements of the solid phase fsiL [8^. l83h . The main difficulties lie in mea- 
suring small signals at low temperatures in the presence oflarge backgrounds. With improv- 
ing experiments measurements were conducted down to 20 inK. While there is still no clear 
evidence of a phase transition in the melting curve experiments, recent pressure measure- 
ments and specific heat measurements have both shown deviations from the expected pure 
Debye lattice behavior. Early on we proposed that these deviations might be related to a 
glass transition and be described by the contributions of two-level systems (TLS). f46Ll48ll 

We model the system of noninteracting TLS with a compact distribution of two-level 
excitation spacings, see Fig. [2] Our results show that the low-temperature deviations in the 
measured specific heat can be explained by contributions from a glassy fraction and/or TLS 
of the solid. 
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3.1 Two level system model for the specific heat 

To avoid complications due to the presence of 3 He atoms, we will compare the effect of 
different growth processes on ultrapure 4 He containing at most (nominally) 1 ppb of 3 He 
atoms. At such low levels of impurities, we expect to see the intrinsic properties of solid. 

We postulate that at temperature much below the lattice Debye temperature, the specific 
heat of solid 4 He is described by 

C(T)=C L (T)+C g (T), (3) 

where the lattice contribution to the molar specific heat is given by C L (T) = B L T 3 , with 
coefficient B L = 12^ 4 /?/5©^, R = 8.314 J/(mol K) is the gas constant, and Q D is the Debye 
temperature. The second term describes the glass contribution due to the TLS subsystem and 
is given by C g (T) = ksR-^f Jq dE @ g (E) f(E), with kg the Boltzmann constant and f(E) 
the Fermi function. The density of states (DOS) of the TLS may be modeled by the box 
distribution function (Fig.[2]i: 

%(E) = ^%[l-tenh((E-E c )/W)}. (4) 

Here S>o is the zero-energy DOS, E c is a characteristic cutoff energy, and W is the width 
of the truncated density of states. For E c — > °°, one obtains the standard hallmark result of 
glasses at low temperatures: C g (T) = B g T, where B g = kgR^o- As we will elaborate in the 
next section, the glass coefficient B g has an intrinsic finite value at low temperature even for 
the purest 4 He samples, independent of 3 He concentration. 

As shown in {4), our model goes beyond the standard glass model by introducing a 
cutoff in the DOS of the TLS (Fig. [2j». The cutoff could be due to the finite barrier height of 
double-well potentials giving rise to the TLS because in real materials the tunneling barrier 
has an upper bound set by lattice and dislocation configurations li^l . At high temperatures, 
the TLS contribution is less important since the thermal energy easily overcomes the barrier 
and effectively resembles a single harmonic degree of freedom. 

3.1,1 Specific Heat 

We compare our calculated specific heat with the experimental data by the Penn State group 
HIH for four different growth processes: BC20, BC04, SL34 and SL31. BC20 (04) is the 
sample grown by blocked capillary (BC) method over 20 (4) hours. SL 34 (31) represents 
the samples in solid-liquid coexistence state with 34 (31) percents of solid ratio. Notice that 
sample SL34 actually corresponds to their reported 75% solid-liquid coexistence sample and 
SL31 corresponds to their constant pressure sample (CP) f49L r77h . The data are described 
with three parameters: 31q, E c and the Debye temperature ®o- We first determine ®o, or 
the lattice contribution, from the high-temperature data. The lattice contribution is then sub- 
tracted from C to obtain the difference 8C = C — Cl. We fit 8C/T with our specific heat 
formula for TLS. 

Next we plot 8C/T in Fig. [3] The TLS model with cutoff describes well the data. In these 
plots we fixed the width of the cutoff to W = 1 jUeV since there is no qualitative difference 
when varying W within reasonable range. Notice that the shape of SC(T) depends strongly 
on the subtraction of the high-temperature lattice contribution. 

The TLS behavior is mainly characterized by the zero-energy DOS and the cutoff energy 
of the TLS, which are both noticeably larger in BC04, see Table l3.lTl This may be explained 
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Fig. 3 SC/T for experiments 
(squares) 17711 and the modified 
glass model with a cutoff en- 
ergy in the TLS DOS (blue line) 
of four samples with different 
structural quality where SC = 
C — Cl- Note the large deviation 
of data points at high tempera- 
tures in the highest purity crys- 
tal SL34 of solid-liquid coexis- 
tence, which makes the subtrac- 
tion of the Debye contribution in 
this sample questionable. 
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by the rapid growth process of a strained crystal, which gives rise to both a larger TLS 
concentration, «tls, and a smaller cutoff energy, E c , i.e., a smaller maximum tunneling 
barrier height. Since the TLS concentrations of these samples range from 3.7 to 21.5 ppm, 
which are at least 1000 times larger than the nominal 3 He concentration, we believe that 3 He 
has no effect on the observed intrinsic heat capacity of ultrapure solid 4 He. 
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24.5 


2.2 


1.7 


3.7 


21.3 


SL31 


25 


21.25 


24.8 


2.9 


2.2 


6.4 


36.9 


BC20 


33 


20.46 


29.7 


3.0 


2.3 


6.9 


39.5 


BC04 


33 


20.46 


28.9 


6.5 


3.3 


21.5 


115.0 



Table 1 Physical and model parameters: Debye temperature ©o, zero energy TLS DOS S>o, cutoff energy 
E c , concentration of TLS «tls = S*o x E c and excess entropy AS. 
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3.1.2 Entropy Analysis 

Our analysis of the excess entropy supports the existence of a glassy component or TLS. 
The excess entropy, 



is associated with an excess specific heat. We find consistently for specific heat measure- 
ments f76Ll55ll77ll that the obtained entropy values AS ~ 10~ 4 J/(K mol) are 5 to 6 orders of 
magnitude smaller compared to the theoretical prediction for a homogeneous supersolid if 
the entire sample actually underwent Bose-Einstein condensation (BEC). In the limit of 
a non-interacting BEC bosons with a quadratic dispersion one finds the standard result 
ASbec = 15/4(£(5/2)/£(3/2))/?(r/r c ) 3 /2 ~ (5/4)7? - 10.4 J/(K mol). This means that 
if AS is indeed due to supersolidity, then the supersolid volume fraction is at most 1 1 ppm or 
0.001 1% in the most disordered or quenched sample of the four ultrapure samples studied in 
this work, i.e., sample BC04. Such a supersolid fraction in the specific heat is more than 100 
to 1000 times smaller than is usually reported for the non-classical rotational inertia fraction 
(NCRIF) in TO experiments. This enormous discrepancy between supersolid fractions in 
specific heat and TO experiments was already noticed in Refs. j4^ . l48ll . while Lin et al. Wfa 
keep invoking a hyperscaling mechanism of unknown origin. Until to date, this discrepancy 
remains a major puzzle that is hard to reconcile within a supersolid scenario. 

The validity of the analysis of the entropy in terms of a non-interacting BEC is repeat- 
edly questioned on grounds of how robust it is in the presence of interactions. We discussed 
in our original study 14611 the effects of interactions on the entropy. Here it is important to 
realize that the entropy is a total count of all low-energy states irrespective of a particular 
model for the specific heat. Hence, we concluded that strong-coupling effects cannot change 
the order of the effect. They may only change the magnitude. For example, the well-known 
strongly correlated superfluid system 4 He possesses a superfluid entropy of ~ 0.67? ~ 4.6 
J/(K mol) I87I1 . which is only half the value of the non-interacting BEC. With hindsight this 
justifies the neglect of strong-coupling effects in the order of magnitude analysis. Clearly 
the entropy is a reliable measure of any phase transition in 4 He. No matter how one looks 
at this puzzle, the reported excess entropy is either far too small to explain observed NCRIF 
effects in torsional oscillators or far too large to describe the boiling off of 3 He atoms from 
dislocation lines, when the nominal concentration of 3 He is less than 1 ppb. For those rea- 
sons, we argued in favor of two-level systems of low-lying states until a better microscopic 
understanding of solid 4 He emerges. 

3.1.3 Comparison with Pressure Measurements 

Next we relate the pressure measurements with the specific heat measurement through our 
model. The quantities to characterize the pressure measurement in the combined lattice and 
glass models are ai and a g defined by 




(5) 



P(T)=P +P L (T)+P g (T)=P Q +a L T 4 +a g T z , 



(6) 



where P(T) is the pressure at temperature T. Pq, Pl, P g are the corresponding pressure 
contributions of the ions at zero temperature, lattice vibrations, and two-level excitations of 



Fig. 4 (P-P ) /T 2 vs. T 2 for P ~ . 
bar. The squares represent the data 1 
ported by Grigor'ev et al. (8ll |8~I 
The blue-thick and black-thin lines a 
predictions for BC04 and BC20, 1 
spectively. The black-dotted line h 
been shifted vertically by a consta 
compared to BC04 to illustrate the c 
pability of describing Grigor'ev's di.._ 
by the TLS model. Predicted curves 
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comparison. 
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The large scatter in data at lowest T 
is due to the subtraction of Pq. Dashed 
lines are guides to the eyes. 
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the glass. On the other hand, the thermodynamic Maxwell relations between pressure and 
specific heat give 

/ v m m 

where Yi are the Griineisen coefficients of the glass (g) and lattice (L). Literature values for 
the Griineisen coefficient of phonons in solid hep 4 He range between 2.6 <Jl< 3.0 f82ll88ll . 
while nothing is known about y g of glassy 4 He. For simplicity we assume y g ~ % = 2.6. 
Equation ^} is related to the first Ehrenfest relation involving the compressibility, which 
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was shown to be always satisfied in glasses l89ll . In Figs. [4] and [5] we show the temperature 
dependence of the (P — Po)/T 2 data reported by Grigor'ev et al. l8lLl82ll and Yin et al. lS7tl . 
The curves for samples BC04 and BC20 are derived from our specific heat analysis of the 
data by Lin et al. 17711 The key result is the dependence ~ T 2 with finite intercept at T = 0. 
In the TLS model, the finite intercept describes the glassy contribution, whereas the T 2 
behavior is attributed to phonons. In conclusion, the data by Grigor'ev et al. and Yin et al. 
are in agreement with predicted (P — Pq)/T 2 curves for a system of TLS. 



4 General response function formalism and physical origin of dynamical anomalies 

To rationalize the TO experiments, we developed a general phenomenological formalism in 
Ref. We have since extended and applied it to visco-elastic and dielectric properties. 
With simple modifications, this predictive approach can be used to study all measurable 
dynamical response functions. Here, we summarize the essence of our approach. In later 
sections, we will apply it in a self-contained way to the various quantities that we wish to 
interrogate. We start with the equation of motion for generalized coordinates g,-. These co- 
ordinates may be an angle, q = 9, in the case of the TO experiments, Cartesian components 
of atomic displacements , qi = m, with u the local atomic displacement (as in our analysis 
of the visco-elastic and dielectric response functions), or any other. Associated with these 
coordinates are conjugate generalized momenta p, (e.g., angular momentum in the TO anal- 
ysis, linear momentum for Cartesian coordinates) and their associated generalized forces Fj 
(torques in the case of the TO, and rectilinear forces for atomic displacements). Physically, 
these forces correspond to a sum of two contributions: (i) Direct forces (Fd, recf ). These may 
originate from either externally applied forces on the solid (F ext ) or lowest order "direct" 
internal forces (fj,»). By direct forces we allude to forces on the coordinate q that do not 
involve the response of the system on q as a result of its change. 

(ii) Indirect "backaction" forces (Fba)- These allude to higher order effects wherein a 
variation in the coordinate q can lead to displacements in other parts of the medium (e.g., 
those involving plastic regions or nearby atoms) which then act back on the original coordi- 
nate q. To linear order 

F BA (t) = f_j(t-t')q(t')dt'. (8) 

The backaction function g(t — t') captures how a displacement q at time t' can lead to a 
perturbation in the solid which then acts back on the coordinate q at time t. With these 
preliminaries in tow, we now outline our standard linear response formalism which further 
takes into account all higher order backaction effects. 

• (1) Write down the Newtonian equation(s) of motion for the generalized coordinate(s) 
q that we wish to study. With Xo a suitable differential operator Xq 1 , these can be cast as 

Xo 1 q(t ) = Fdirect (t ) + F BA (f ) . (9) 

This equation might seem a bit formal. To make it concrete, we note that in the simplest case 
that we will discuss- that of the compact scalar angular coordinate 9 for the TO orientation, 
the operator 

_i d 2 d 

*0 = r osc -7-T + Jose -77 + dose (10) 
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where I sc,Yosc, and a osc are, respectively, the oscillator moment of inertia, dissipation, and 
torsion rod stiffness. A more complicated tensorial operator involving the elastic modulus 
appears when writing the equations of motion for the Cartesian components of the atomic 
displacements. 

• (2) When Fourier transformed to the complex frequency (co) plane, Xo{g>) corresponds 
to the bare susceptibility. We will denote the Fourier transforms of the various quantities 
(forces, backaction function) by simply making it clear that the argument of the various 
quantities is now the frequency co and not the time t. We trivially recast Eq. ^ as 

X~ 1 {C0)q(C0) = F direc ,(C0), (11) 

where 

X'Hco)=X 1 (co)-g(co). (12) 

Equation d 1 2b will be used in all of our upcoming analysis. The physical content of Eq. 
can be seen by writing its inverse as a geometric (or Dyson) series 

X(t)=Xo(t) + fdt'xo(t)g(t-t')Xo(t') 

+ jdt'J df' X o(t)g(t-t') X o(t'W-t")Xo(t") + ~~ (13) 

The terms in this series correspond to (a) the direct contribution (Xo), (b) a lowest order 
backaction effect wherein a displacement at an earlier time t' <t leads to a deformation of 
the solid which then acts back on the coordinate at time t (the term J dt'xo (t)g{t — t')Xo (t')), 
(c) a higher order process in which a deformation at time t" leads to a backaction from the 
solid on the coordinate at time ?' > t" which in turn then acts back on its surroundings 
which then acts back on the original coordinate at time t (the term / dt' J dt"Xo{t)g(t — 
t')Xo{t')g{t' — t")Xo(t")), and so on ad infinitum. In Fourier space, the convolution integrals 
become products and Eq. d 1 3b becomes 

Z(ffl) = Xo(co) +Xo(a)g(®)Xo(®)+Xo(a>)g(co)Xo(co)g(co)Xo(co) + (14) 

The geometric series of Eq. d 1 4b sums to Eq. d 1 2b - As simple as it is, Eq. d 12b combining 
standard linear response theory (step 1) with backaction effects (or the Dyson equation) of 
step 2, leads to a very powerful tool that allows us to investigate numerous systems while ac- 
counting for arbitrarily high order backaction effects. As is well known and we will expand 
on and employ in later sections, the real and imaginary parts of the poles of the suscepti- 
bility x( tt> ) allow us to probe typical oscillation times and dissipation. This will allow us to 
connect with TO and other measurements and make precise statements about the backaction 
function g which affords information about the dynamics within the solid. It is important 
to stress that Eq. ( 112b is general. In the derivation above no assumptions need to be made 
concerning the precise physical origin of the backaction function g. The adduced function 
g captures all effects not present in the direct equations of motion for the normal solid. If 
supersolid effects would be present, they will directly appear in the function g. 

Thus far, our sole assumption was that the deformations q are small enough to justify the 
linear (in q) order analysis of Equations ([8]|9]l for measurements on the rigid solid. We now 
invoke additional assumptions (which have been partially vindicated in a growing number of 
experiments since our original proposal j47ll ). These assumption relate to the form of g(co) 
and its dependence on temperature. They are motivated by our view of defect quenching and 
characteristic relaxation times as the origin of the 4 He anomalies. 
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• (3) As we will elaborate on in later sections, data for disparate susceptibilites % taken 
at different frequencies CO or temperatures T collapse onto one curve. This indicates that g 
is a function of only one dimensionless argument (C0T(r)) instead of both CO and T inde- 
pendently. That is, there is only one dominant temperature dependent relaxation time scale 
t(T) for the backaction of the quenched solid on the original coordinate q. As is well known 
and alluded to in Section 10, exponential damping with a single relaxation time T, leads to a 
function g S mgie{ t0 ) = £o(l — )/(l + C0 2 T 2 ) which, when plotted with the real and imag- 
inary parts of g along the horizontal and vertical axes describes a semi-circle as a function 
of the dimensionless quantity (cot) > 0. However, when plotted in this way, the 4 He data 
for the complex susceptibility measured by TO and other probes lead to a collapsed curve 
which is more like that of a skewed semi-circle and there is a distribution of relaxation times 
about a characteristic time scale T. For ease of analysis, we will approximate the complex 
response functions by Eqs. (0. The curve collapse allows for information about how the 
characteristic transient relaxation times t(T) increase as T is decreased. We further invoke 
the Vogel-Fulcher-Tammann form for glasses l9(ih . 

f Tne A /( T ~ T °} for T > 7n 

Zt<1: (15) 

Here, 7b is the temperature at which the relaxation times would truly diverge and A is an 
energy scale. In fitting the data in this way, negative values 7b < were often found. That 
is, the typical dynamics as adduced from the collapsed t(T) is faster than simple activated 
dynamics (one in which 7b = in Eq. d!5ll). This is consistent with the intrinsic quantum 
character of the solid 4 He crystal with large zero point motion as compared to classical 
activated dynamics. 

What is the physical content of this general formalism vis a vis the putative supersolid 
transtition? Given perturbations of a typical frequency CO, the backaction response g from the 
plastic components acting on q may either be sufficiently rapid or slow to respond. Just at the 
tipping point when cot(Tx) ~ 1 different components will be maximally out of synchrony 
with each other in being able to respond to the perturbation or not and the dissipation (given 
by the reciprocal of the imaginary part of the zero of % _1 (co)) is maximal. Similarly, there 
will be a change in the typical periods of the system between a system which at high T (i.e., 
T < Tx) contains rapidly equilibrating plastic components to those at low T, which are too 
slow to respond and thus the system appears to have undergone "a transition". 

In the sections that follow, we will apply and replicate anew and at great length the con- 
siderations outlined above to the particular set of physical quantities that we wish to investi- 
gate. We start, in Section[5]with the investigation of the TO (for which the above formalism 
was first developed) and then move to explore other arenas- the viscoelastic (Section [6} and 
dielectric response (Section0 functions where the above formalism leads to experimentally 
testable predictions. 



5 General susceptibility and response function of torsional oscillators 

A formulation of the rotational susceptibility of the TO was given in Ref. |47h . It is now often 
used as a basis for the linear response discussion l9ll[l4l[T5l . In Section © we outlined the 
key points of this formulation when written in its general form. Since its derivation it has 
been applied by us and others to study other response functions. The result of Eq. i ll St below 
is none other than Eq. d 1 2b when the generalized coordinate q corresponds to the TO angle. 
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The bulk of this section will be devoted to analyzing the experimental consequences of this 
relation and its related counterpart for the double TO. 

It is important to realize that the TO experiment measures the period and dissipation 
of the entire apparatus by reporting the relationship between the force and displacement 
(angle). Therefore a model is needed to determine the relation between observable period 
and dissipation and the moment of inertia, damping and effective stiffness of the media. We 
start with the general equation of motion for a harmonic TO defined by an angular coordinate 
6 in the presence of an external and internal torque. l47ll 

I osc 0(t) + y scd{t) + a osc e(t) = M ex ,{t) +M int {t). (16) 

Here, I osc is the moment of inertia of the (empty) TO apparatus, a osc is the restoring constant 
of the torsion rod, and y osc is its damping coefficient. M ext (t) is the externally imposed torque 
by the drive. A/,-„,(f) = / g(t — t')8(t')dt' is the internal torque exerted by solid 4 He on the 
oscillator for a system with time translation invariance. In general, the backaction g(t — t') 
is temperature, T, dependent. 

The external torque, M ext (t) = L(t), is the derivative of the total angular momentum of 
a rigid body, L(t) = 4/ cPx p(x)r 2 6(x), where r is the distance to the axis of rotation, 
p(x) is the mass density and 0(x) the local angular velocity about the axis of rotation. I92I1 
The experimentally measured quantity is the angular motion of the TO - not that of bulk 
helium, which is enclosed in it. Ab initio, we cannot assume that the medium moves as 
one rigid body. If the non-solid subsystem "freezes" into a glass, the medium will move 
with greater uniformity and speed. This leads to an effect similar to that of the nonclassical 
rotational moment of inertia, although its physical origin is completely different. We argue 
for an alternate physical picture, namely that of softening of the oscillator's stiffness. 

The angular coordinate 9(t) of the oscillator is a convolution of the applied exter- 
nal torque M ext {t) with the response function %{t,t'). Causality demands = 6(t — 
t')x(t,t')- Under Fourier transformation, this leads to the Kramers-Kronig relations. In any 
time translationally invariant system, the Fourier amplitude of the angular response of the 
TO is 

X l (co)9(co)=M ext (co)+M inl (co). (17) 

Defining the total angular susceptibility as x~ l = Xq 1 —Mi„,, we write the effective inverse 
susceptibility as 

X~ l {®) = a osc -iY osc co-I osc a 2 -g{co), (18) 

where g(co) is the Fourier transform of the backaction due to the added solid 4 He. In what 
follows, we will treat the backaction as a small perturbation to the TO chassis. 

We will now apply our formalism to the study of the ingle TO, which is described 
by Eq. d 1 8b . and then turn to the double TO. Very recently, Beamish f93tl and Maris l94ll 
employed the same general linear response formalism to explain some of the TO results in 
terms of purely mechanical effects due to either the changing stiffness of the torsion rod or 
Soppiness of the sample cell flange (lid). 

5.1 A model for the single torsional oscillator 

In what follows, we analyze the experimental consequences of Eq. 1 II8I ). 
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5.1.1 Rotational susceptibility - period and dissipation 

We can now calculate specific consequences of the phenomenological model introduced 
above. The effective oscillator parameters are defined as the sum of parameters describing 
the apparatus, Xq < an d tne added solid 4 He given by 

g(ffl) = iYHe(0+lHe(0 2 +gaG{co). (19) 

It is convenient to introduce a net moment of inertia / = I osc + In e and net dissipation 
7= Yosc + Yhc- The transient dynamic response function G(co) can be approximated by a 
distribution of overdamped oscillators with relaxation time T as discussed in Section $2$ 
[see Eqs. ©, in particular]. 

The resonant frequency of the TO with a backaction is given by the root of 

X~ l {<o) = a-iyco-Ico 2 -g G(co) = 0. (20) 

As discussed in Section [4] we anticipate that when the relaxation time is similar to the 
period of the underdamped oscillator, the dissipation will be maximal, sometimes referred 
to as "cot = 1" physics. 

In linear response theory, the homogeneous Eq. d20b is scale invariant. Thus, we normal- 
ize all oscillator quantities by the effective moment of inertia /, i.e., a = a/I, y = y/I, and 
go = go/ '!• As can be seen from Eq. i !20b . for an ideal dissipationless oscillator (y = 0), the 
resonant frequency C0q = \/a~ is the pole of %{co) in the limit 1/t — > 0. If we expand % 
about this root, CO = COq + 8 CO, we find to leading order in 8 CO 

if +2(00 

It follows that the shift in dissipation with respect to high temperatures is 

AQ-^Q-'-Q^^^ImGicoo), (22) 
co 

whereas the shift in resonant frequency is 

Aco = 2n{fo-f)^ -^(2c*oReG(c*o) + 7lmG(c*o)), (23) 

which increases monotonically when T is lowered. Combining Eqns. ( 122b and d23b for the 
strongly underdamped oscillator, we arrive at 

AQ 1 _ 4lmG(coo) _ 2 Im G(co ) 

Aco 2c*oRe G(coq) + ylm G(coq) ~ coo Re G(coq) ' 

It is this general relationship for the response function of the damped oscillator that was 
successfully applied in the TO analysis by Pratt et al. fflfisll to demonstrate the interplay 
of rotational, relaxation, and shear dynamics in solid 4 He. For a Debye relaxor the ratio of 
Eq. 024b reduces to 2t and provides a direct means to measure the relaxation time. Similar 
results for the ratio were obtained for other phenomenological models. I65l . l9lll For example, 
Huse and Khandker l65ll assumed a simple phenomenological two-fluid model, where the 
supersolid is dissipatively coupled to a normal solid resulting in a ratio of AQ~ l ^ « 1, 
while Yoo and Dorsev l9lll developed a viscoelastic model, and Korshunov|39] derived a 
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Fig. 6 The resonant frequency (black, left axis) and dissipation (red, right axis) vs. temperature. The experi- 
mental data from Hunt et al. flHl are well described by a Cole-Cole (CC) distribution function. I95ll96ll 




Fig. 7 Cartoon of the double torsional oscillator modeled in Eqns. )25t . The upper moment of inertia (I[ ) is 
the dummy bob, while the lower moment of inertia (Ii) is the cylindrical pressure cell that can be loaded with 
4 He. The stiffness of the torsion rods is given by k[ and ki with k[ £2 by design. 



two-level system glass model that captures the results of the general model originally pro- 
posed by Nussinov et al. f47h To make further progress we assume that T follows the phe- 
nomenological Vogel-Fulcher-Tammann (VFT) equation of Eg. d 1 5 b - 

Figure|6]provides a fit to the measured data by Hunt et al. 113(1 assuming a CC distribution 
of relaxation times. As shown, an excellent fit is obtained. For comparison, we also tried a 
DC distribution for relaxation times, but found only fair agreement. It is worth mentioning 
that unlike in the Debye relaxation analysis by Hunt and coworkers, i.e., a single overdamped 
mode, we do not require a su pers olid component to simultaneously account for frequency 
shift and dissipation peak l95ll96ll . Our model leads to a universal scaling of period change 
vs. dissipation in a Cole-Cole or Davidson-Cole plot as seen in Ref. fl^LITsll. Indeed similar 
viscoelastic behavior may have already been observed in solid hydrogen 19711 . 



5.2 A model for the double oscillator 



The double TO results of the Rutgers group by Kojima have proved difficult to explain when 
simply extrapolating from the single TO model [91]. Here we model the coupled double TO, 
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T(mK) 

Fig. 8 Frequency and dissipation in double TO by Aoki et al. floll (symbols) compared with glass theory 
(lines). Panel (a): Temperature dependence of resonant frequency shifts Af\ (black, left axis) and Afi (blue, 
right axis). Panel (b): Temperature dependence of dissipation Q[ 1 (black, left axis) and 1 (blue, right axis). 
The experimental data were corrected for the significant temperature-dependent background of the empty cell. 

sketched in Fig. [7J by the following system of equations: 

(-ha 2 -iy x w + k x + k 2 ) © 1 {w)-k 2 @2{(o) = F(a>), 
(-I 2 a 2 -iY20) + k2-g{co))©2(ca)-k2©i{(o) = 0. (25) 

where ©, are torsion angles, Yi are damping coefficients, are torsion rod stiffnesses, g((o) is 
the glass backaction term, and F(co) is the applied external torque. The subindex "1" refers 
to the upper or dummy bob in the experiment, while "2" refers to the lower oscillator with the 
pressure cell that can be loaded with solid 4 He. For a strongly underdamped oscillator and a 
small backaction, it suffices to solve first for the bare resonant frequencies and later include 
perturbatively damping and backaction terms, for details see Ref. l95ll . More recently this 
approach has been extended to a triple TO f9^1 . 

Figure [8] shows good agreement between our phenomenological model of the coupled 
double TO and experiment. The TO parameters /, and k{ can be determined from the bare 
resonant frequencies jf = (of /2k. In addition, the damping coefficients ji can be extracted 
from the high-temperature dissipation Finally, the backaction g(co) accounts through 
t(T) for the temperature dependence of Af, and Q[ . Our phenomenological theory of the 
double oscillator explained for the first time both frequency shift and dissipation peak for 
in-phase and out-of-phase torsional response fioh . Data for in-phase frequency f\ = 496 Hz 
and out-of-phase f% = 1 173 Hz are shown in Fig. [8] plotted against the temperature. The 
obtained values for moment of inertia and rod stiffness agree well with other estimates |9^| . 

Remarkably an anomalous damping coefficient /i ~ — 72 is required to explain the be- 
havior of increased dissipation with increased frequency. Such anomalous damping is al- 
ready required to describe the unloaded pressure cell. Thus it is unrelated to the properties of 
solid 4 He and an intrinsic property of the composite TO. After loading the cell with solid 4 He 
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the dissipation ratio becomes £?2 /fii = 2-5 at 300 mK with frequency ratio /2//1 = 2.37. 
Our fit results in a negative parameter 7o = —32.73 mK. This effective negative value of To 
is in line with earlier comment in Section[4] concerning the quantum character of solid 4 He. 
This value may be indicative of strong zero point quantum fluctuations that thwart a glass 
transition. 

Finally, the comparison in Fig.[8]shows that an explicit frequency-dependent backaction 
must be used with go(co) = go (a)/(of] and p = 1.77 to account for the experimental fact 
of Af\jf\ « Af%lf%, i.e., the relative frequency shift is unaffected by the changing reso- 
nant frequency. In contrast, various theories describing solid 4 He in torsional oscillators as 
viscoelastic material l9lll or two-level systems moving through a solid matrix I37L l38l . l39ll 
predict an exponent of p = 4 for the backaction term. 



6 Shear and stiffness of a viscoelastic medium 

Another aspect of the d y nam ic respons e of 4 He crystals was revealed through a series of 
elasticity studies . I22I. 1201 1 1 Od , flOll [5^1 In particular, Beamish and coworkers demon- 
strated the qualitative similarity between shear modulus and the TO measurements. In the 
shear modulus experiment solid helium is grown in between and around two closely spaced 
sound transducers. When one of the transducers applies an external strain, the other trans- 
ducer measures the induced stress from which the shear modulus of the sample is deduced. I lOOl l 
In this way, the experiment provides a direct measurement of the elastic response to the ap- 
plied force within a broad and tunable frequency range. The frequency dependence is espe- 
cially crucial in determining the nature of the relaxation processes and complements current 
TO experiments with their limited frequency range. 

Similar to the TO analysis in the previous section, we analyzed the shear modulus within 
the general linear response function framework, where the amplitude of the shear modulus 
increases (stiffens) upon lowering T, because of the freezing out of low-energy excitations. 
This change is accompanied by a prominent dissipation peak, indicative of anelastic relax- 
ation processes. We calculated the complex shear modulus [i(co; T) of a viscoelastic material 
and predicted: (a) the maximum of the shear modulus change and the height of the dissipa- 
tion peak are independent of frequency and (b) the inverse crossover temperature l/Tx vs. 
the applied frequency © obeys the form 0)t(7x) = 1 characteristic of dynamic crossover. 



6. 1 Model of dynamic shear modulus 

As in our analysis of the TO, we start with the same general linear response function formu- 
lation outlined in Section[4] Our final result of Eq. J29l l will, once again, reflect the general 
relation of Eq. d 1 2b - Here we replace the angular coordinate of the TO with a displacement 
coordinate and the restoring force with a stress tensor l50ll . For ease, our notation in the 
below will differ slightly from that in Section^ The equation of motion for displacement m, 
in the ;-th direction of a volume element in the presence of an external driving force is 

-patut+djojj* =f? XT (co)+f? A (co), (26) 

where p is the mass density, ff XT and /) BA are the external force density and the backaction 
force density, and <jfj e is the elastic stress tensor of solid helium. In general, ofj e = Aya uy, 
with the elastic modulus tensor A,-^/ Il02ll . In the case of a homogeneous solid with shear 
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wave propagation along the z axis and wave polarization in the x-y plane, the backaction 
takes on the form 

f? A (w) = G(w,T)d? Ui (w), (27) 

where G describes the strength of the backaction on the solid (viscoelastic response) and 
i = x,y. Although fP A (co) is much smaller than the purely elastic restoring force dj ajj e , it is 
this term that is responsible for the stiffening of shear modulus with decreasing temperature. 

Polycrystalline and amorphous materials are nearly isotropic, hence the elastic modulus 
tensor becomes Ayj; = AoS/jS^ + ^0(8^811 + 5//5^). Note the stress tensor in Eq. ( 126b is 
finite only for orientations j = z and either k or / equal to z- With k, I being interchange- 
able, the relevant element will be A/ z ,- z , which gives the purely elastic shear modulus p^. 
Finally the fully dressed shear modulus (dressed by the backaction) relates the displacement 
to an external force, or [— p CO 2 + p.d 2 ]uj(co) = /^ XT (co). Comparing this expression with 
Eq. ( 126b . we obtain for the dynamic shear modulus in a viscoelastic material the general 
response function 

t i{<o;T) = iM ) {T)-G{(d;T). (28) 

Next we employ for G(co; T) the Cole-Cole distribution function [specifically, by reference 
to Eq. ©, we set G(co; T) = gofioG(co, T) to obtain the specific form 



pi{co;T) = j&o 



1 80 



1 - (/COT) 



(29) 



with the sample dependent parameter go- The experimentally measurable quantities are the 
amplitude of the shear modulus, |ju|, and the phase delay between the input and read-out 
signal, (j> = arg (/x); (j) measures the dissipation of the system, which is related to the inverse 
of the quality factor g 1 = tan0. 

Several interesting results follow immediately from the general response function in 
Eq. J29b : (1) The change in shear modulus A p. between zero and infinite relaxation time 
is Afi/pio = go- It measures the strength of the backaction as well as the concentration of 
defects. (2) At fixed T, the shear modulus amplitude |ju(co; T) \ decreases with increasing CO. 
(3) The parameter cot is the only scaling parameter. (4) The peak height A(j) of the phase 
angle is proportional to go- When cot ~ 1, then A<$> « gocot(o!7r/4)/(2 — go) for go <S 1. In 
the limit 1 < a < 2, this simplifies further to 

40«(l-a/2)(4/i//io)<l, (30) 

where A (j) is in units of radians. Quite remarkably, the peak height A(j) depends only on the 
phenomenological parameters a and g °c A /A. At fixed temperature the maximum change of 
both Afl and A<j) is independent of frequency. 



6.2 Results 

Let us compare our model calculations with the experimental shear modulus measurements 
by Beamish and coworkers B^l for a transducer driven at 2000 Hz, 200 Hz, 20 Hz, 2 Hz, 
and 0.5 Hz jl03ll . Specifically, for the T-dependent t(T) in Eq. d29b we consider Vogel- 
Fulcher-Tammann (VFT) and power-law (PL) relaxation processes. In our model parameter 
search, we are not constraining 7q to be positive (where T diverges). Fair agreement between 
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Fig. 9 Experimental data and theoretical calculations of the shear modulus vs. temperature assuming a VFT 
relaxation time. The red and blue squares are the experimental data for the amplitude and dissipation of shear 
modulus. The black-solid lines show the theoretical calculations. We used the set of parameters a = 1.31, 
g = 1.44 x 10 _1 ,A) = 0.47 pA Hz" 1 , T = 50.0 ns, A = 1.92K, and 7b = -69.3 mK. Notice a negative 
Tq means that there is no true phase transition occurring at finite temperatures, probably because of strong 
quantum fluctuations of helium atoms. The shear modulus amplitude and the phase angle were shifted by 0. 1 
pA Hz - 1 and 3 degrees with respect to the 2 Hz data. 



calculations with a single set of parameters and experiments is obtained with 7b = —69.3 
mK, see Fig. [9] We refer to these calculations as "VFT<". Similar to the TO results, a neg- 
ative 7b means that no real phase transition occurs in solid 4 He. The agreement between 
theory and experiment at highest (2000 Hz) and lowest (0.5 Hz) frequencies is of poorer 
quality. This discrepancy may be related to very different backgrounds at these frequencies 
and may not be intrinsic to 4 He. Another reason may be the presence of additional dissipa- 
tion mechanisms, not accounted for. Furthermore, our calculations confirm that the change 
of amplitude and peak height of phase angle are nearly frequency independent between 2 
Hz and 200 Hz. By defining the crossover temperature Tx as the point where the phase angle 
peaks, we find as predicted that Tx decreases with decreasing ft). 

Finally, when we set 7b = K the VFT expression ("VFTo") reduces to an activated 
Arrhenius form. While it describes reasonably well the position Tx, it gives a much narrower 
linewidth for A(j) than does VFT<, which is not in accord with the data and thus not shown. 
Notice that the VFT relaxation is not the only possible relaxation process that can describe 
the data; power-law or other types of relaxation can give similar level of agreement to the 
experiment I56I1 . Iwasa proposed a relaxation process Il04ll similar to our phenomenological 
one, which is based on the theory of pinned vibrating dislocation lines by Granato and Liicke 
Il05ll . Clearly further experiments at lower frequencies and lower temperatures are required 
to determine the exact type of relaxation processes in solid 4 He. 
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Fig. 10 The Cole-Cole plots for exper- 
imental data and for VFT calculation. 
For given form of T, all different fre- 
quency curves collapse onto one single 
master curve reflecting that COT is the only 
scaling parameter. The Cole-Cole plots 
show reflection symmetry about Re[|/i — 
fk)\/A tl]=0.5, which is a consequence of 
the Cole-Cole distribution function. 
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Fig. 11 Prediction for the inverse crossover 
temperature vs. applied frequency. The 
green squares correspond to Tx in Ref. ($(\ \. 
For the power-law process with phase tran- 
sition occurring at 40 mK, we used T = 
M\T g \/(T - To))" for T > T and T = ~ 
for T < T . 
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Figure [TO] shows the Cole-Cole plot for experiments and calculations. The experimental 
curves for different frequencies collapse roughly onto one curve confirming our theoreti- 
cal assumption that cot(T) is a universal scaling parameter. This behavior was also seen in 
TO experiments fl3h . In addition, the data are symmetric about Re[|,u — Aji] = 0.5 as 
expected for the Cole-Cole distribution. A more detailed data analysis may resolve the re- 
maining discrepancy between theory and experiment. The discrepancy may be due to either 
the presence of additional relaxation processes at temperatures above Tx, i.e., a more com- 
plicated form for l(T), or by a modified functional form of the defect distribution function. 

The pertinent question about a true phase transition at zero frequency vs. crossover 
dynamics can be addressed by investigating the asymptotic limit of cot(T) = 1. From this 
expression we estimate Tx as a function of the applied frequency / = CO /2k. Figure. [TT1 
shows l/Tx vs. /. The VFT< and VFTo calculations give significantly better agreement 
than the PL calculations with phase transitions occurring at K (PLo) and 40 mK (PL40). 
For positive 7b (see PL40), we find a true freeze-out transition, which would indicate arrested 
dynamics for / — > Hz. For both VFT and PL relaxation times our calculations demonstrate 
that in the low-frequency limit the existence of a phase transition shows clear signatures 
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Fig. 12 Lump circuit of the generalized 
Maxwell model. The elastic spring Hq rep- 
resents the purely elastic shear modulus at 
high temperature. Each Debye relaxor is 
made out of a series of elastic spring (Irs 
and dissipative dash-pot r/'"'. 




of Tx converging toward the ideal glass temperature Tq. Therefore the absence of arrested 
behavior may serve as experimental evidence against a true phase transition. 



6.3 Viscoelastic model 

The viscoelastic model successfully describes composite materials with anelastic properties. 
In fact, Yoo and Dorsey £911] applied the viscoelastic model to the TO experiments. More 
generally, a distribution of viscous components, embedded in an otherwise elastic solid, may 
be treated through a generalized Maxwell model 1501 15811 . Conceptually one may subdivide 
the material into many elements and solve the coupled materials equations for stress and 
strain. Here we use constitutive materials equations to show that our results for a glass 
are equivalent to the generalized Maxwell model with parallel connections of an elastic 
component with an infinite set of Debye relaxors, see Fig.[T2l 

The anomalous stiffening of the shear modulus can be described within t he vis coelastic 
model, though other defect mechanisms like dislocation glide are possible too ll06Ll59ll . The 
equivalent lump circuit model is sketched in Fig. [12] where the basic dissipative element is 
the Debye relaxor. It is comprised of a serial connection of a rigid solid (RS), characterized 
by an elastic shear modulus /Irs, and a Newtonian liquid (NL) or dash-pot, characterized 
by a viscosity r\ . The RS component describes the ideal elastic solid helium of this volume 
element, while NL represents the anelastic component, which gives rise to viscous damping. 
The two parts are connected in series, so that both share the same magnitude of stress, while 
the net strain is additive. The strain rate equation for both constituents of the Debye relaxor 
is 

e = <7/^RS + <r/Tj, (31) 

where e is the net strain of the Debye relaxor (DR) and a is the magnitude of stress shared by 
the components RS and NL. In order to obtain the above equation, we used the constitutive 
materials equations for strain and stress: £rs = Ors/^Irs and £nl = OVl/tj. After performing 
the Fourier transformation we obtain for a single Debye relaxor (DR) the shear modulus 
/i D r = a/e, 
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with relaxation time Tdr = T] /Mrs- For a viscoelastic material with a single relaxation time, 
the solid behaves like a parallel connection between the elastic part and the Debye relaxor 



H(a>) = Mo+Mdr(o)) = Ao 



Ho 



1 - iCOTuR 



(33) 



with go = Mrs /Po an d Mo = Mo + Mrs is the dressed elastic shear modulus. To consider the 
general case of many components, we introduce Debye relaxors with different relaxation 
times connected in parallel as shown in Fig. [12] The total anelastic contribution from n 
such constituents is given by a weighted sum. The corresponding continuous version of this 
expression with a distribution P(s) of relaxation times is 



Ma^to) = Mrs / dsP(s) 
Jo 



1 



1 - icoxs 



Mrs 



Mrs 



1 — (icox) c 



(34) 



Here the Cole-Cole distribution, 



P(s) 



-(!-«) 



sincuTT 



\+s 2a + 2s a cos ait ' 



(35) 



was used. Similar to our response functions elsewhere, the net shear modulus of the compos- 
ite system is given by the sum of two terms- the purely elastic component and the anelastic 
contribution 



/i(g>) = M)+M«e(») =Mo 



go 



1 - (icox) 



(36) 



Indeed, this expression is identical to the one obtained previously using the general response 
function formalism with a backaction d 1 2b . This is no coincidence, since the backaction term 
accounts for damping and thus describes the anelastic response of defects to the external 
stress. 



7 Dielectric properties of a viscoelastic medium 

The measurements of e(co,T) by Yin et al. J57ll showed the anomalous increase of the 
dielectric function of solid 4 He at low temperatures. A similar test experiment in liquid 
helium showed no such effect. We propose that these results may be explained by an electro- 
elastic coupling of a quenched solid with frozen-in internal stress. Such behavior cannot 
be described by the standard Clausius-Mossotti equation via a change in mass density or 
polarizability (due to, e.g., dipole-induced dipole interactions). Neither the measured change 
of the mean mass density 8p / p ~ 10~ 6 , nor the predicted co rrection in polarizability, which 
actually leads to a decrease of e(ft), T) at low temperatures Il07tll08h . can account for the 
reported anomalous change of the dielectric function of order Se/e ~ 10~ 5 , while a model 
with frozen-in stress and electro-elastic coupling can explain the data. 
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7. 1 Model for electro-elastic properties 

The relation of Eq. d45b that we derive below constitutes yet another realization of our gen- 
eral relation of Eq. dl2b . We now turn to the specifics of the electro-elastic coupling describ- 
ing the interaction between the electromagnetic and the strain fields. The coupling may be 
obtained by expanding the dipole moment, p(r fl ), of a nonpolar atom around its equilibrium 
value: 

p(r fl )«p(R £ ,) + (u fl -V)p| Ra , (37) 

where r a = R fl + u fl is the position of atom a, R a is its equilibrium position, and u„ is 
its displacement. The polarization is obtained by averaging p over a macroscopic volume 
element v, P(r) = (l/v)£ a p(r a ). In the continuum limit (when the macroscopic length 
scale is far larger than the atomic length) £ fl p(r fl ) w (1/v) f v dr' p(r') . In the presence of a 
local strain field, to linear order in the displacement, Eq. (137b reads 

P(r) « (1/v 2 ) JdR [p(R) + (u- V) p(R)]. (38) 

Integration of the first term yields the macroscopic polarization for zero internal strain (a 
solid in equilibrium), Po. It is related to the macroscopic electric field E by Po = XoE where 
Xo = (£o — 1)/4k is the zero-strain susceptibility and £o is the permittivity. The second 
term in the integration describes the polarization change 5P due to atomic displacements. 
Neglecting surface contributions the second term modifies the polarization P = Pq + 5P by 

5P= -(1/v 2 ) J dR (V-u) p(R) « -e s P , (39) 

with en = ( 1 /v) / dR (V-u) the macroscopic frozen-in dilatory strain (we use the repeated 
indices summation convention eu = en + e22 +£33)- In obtaining Eq. <39t , we assumed 
that P is slowly varying. This long-wave length approximation holds for wavelengths of 
microwave (A > 10~ 3 m) and above. The polarization change alters the dielectric function 
= {eoja + Sea by 

Ses (co, T) = -Anxo e u (a, T) (40) 

with Eu — 1 = AKPi/Ei. Only the diagonal components of the strain tensor play a role in a 
leading order expansion of the electro-elastic coupling. In principle, the shear strain can cou- 
ple to the electric field by considering dipole-induced dipole interactions (van der Waals), 
which is a higher order effect. The same coupling mechanism between photon (electric field) 
and phonons (strain) gives rise to the acousto-optical effect. Our derivation of Eq. d40t is, to 
leading order, equivalent to the Pockels coefficient for acousto-optical coupling in isotropic 
dielectrics, 5e,y = — e 2 [2P44 e u +P12 £kk^ij\, where are the reduced Pockels coefficients 
(of order unity) flQgL [TlcL [TTTn . 

We next turn to the locally frozen-in strains and write the equation of motion within the 
general response function theory. As before in an isotropic medium the elastic stress tensor 
ofj e = Xijki du.fr/dxi with Xijft = XoSijSu + lM){§ik8ji + Sii$jk)- If the electric field couples to 
local density fluctuations only through dilatory strain, then the important matrix element is 
the Lame parameter Ao- We write the displacement to an out-of-equilibrium internal (INT) 
force in the presence of the backaction as 



- P co 2 Ui{co) +Xod? Ui (w) = ft m (to) + f? A (co), 



(41) 



26 



where Uj is a displacement in the zth direction and p is the mass density. The backaction 
force density of nearby atoms is given by 

ff\(o) = G{(0;T)df «,-(©), (42) 

where G is the strength of the pertinent backaction. ff NT is the out-of-equilibrium (frozen- 
in) internal force density in the z'th direction at the defect. Integrating Eq. d4 1 b yields the 
strain due to the internal dilatory stress ct,,-: 

e a {(0,T) = CT,.(A -G(co;r)) _1 . (43) 

Again, we assume that the backaction can be described by a Cole-Cole distribution of Debye 
relaxors, G{(0\T) = goAo/[l — (/cot)"] . The corresponding local dilatory strain reads 

e a {(0,T) =e (l-g /[l-(icoT) a }y\ (44) 

where eo = Ou/^o at J = 0. From Eqns. l !40t and(!44b the change in the dielectric function 
due to local strain fluctuations is 

8e ti = -4Tzxoe Q {l - go / [1 - (icoz)"})- 1 . (45) 

This result is similar to the one for the TO and shear modulus discussed in the earlier sec- 
tions. At low temperatures, % — > °° and e„ — > eo, hence the strain is minimal and the reduc- 
tion of the dielectric function due to local strain fluctuations is small. At high temperatures, 
T — > and e„ — > eo(l — go) reaches its maximum resulting in the largest reduction of 
strain, where solid 4 He is softest. The main result is that the dielectric function reflects the 
arrested dynamics of the glassy components at low temperatures through the acousto-optical 
(electro-elastic) coupling. 

The derivation of Eq. i45 b for the dielectric response is extremely general. We illustrate 
how a Cole-Cole form for the elastic relaxation implies a similar dielectric response (and 
vice versa). An identical result holds for other forms of the local elastic relaxation dynamics- 
these will leave a similar imprint on the dielectric response. Historically (since the 1940s) 
the Cole-Cole response function f7lll was found to be valid for the dielectric response. In our 
initial works on 4 He in trying to capture the local quenched dynamics, we first assumed this 
form for the TO and then for the general elastic response. By virtue of their inter-relation 
and coupling, the local relaxation dynamics is the same for all of these quantities. Therefore 
measurements of the dielectric response may inform about local dynamics and vice versa. 
Practically, our prediction of Eq. d45b applies to any nonpolarized system with a local dis- 
tribution of stress relaxations. In polarized materials, detecting a change in the dielectric 
function due to atoms sensing different local fields is challenging because the large intrinsic 
polarizability of the material will overwhelm the contributions derived above. Among the 
nonpolarized solids, solid 4 He favors the observation of this phenomenon especially because 
of its softness. Rapid cooling of solid helium allows a large local strain build-up, which is 
proportional to the size of the effect. Similarly, we expect the effect to be seen via delicate 
effects in solid 3 He, hydrogen or xenon. 
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Fig. 13 Experimental data and calculations 
of the dielectric function vs. temperature. 
The red circles are the experimental data 
of the dielectric function (data by Yin et 
al. (5lh ). The black lines are the calcu- 
lated amplitude and phase lag (dissipation) 
of e(co;T). We used parameters a = 1.49, 
e = 8.88 x l(T 4 ,g = 0.21, T = 10.4 ns, 
A = 1.92 K, and 7b = -119 mK. 
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7.2 Results 



The results of our electro-elastic predictions for the dielectric function are shown in Fig.[l3l 
We obtain excellent agreement with experiment l57ll for an applied alternating voltage at 
500 Hz. Our analysis predicts, that similar to the TO and shear modulus, a dissipation peak 
appears in the dielectric function. The phase lag angle (j) = arg(e) records the lag between 
the real and imaginary part of e(co;T). Future observation of the dissipation peak will pro- 
vide an important test of our picture concerning quenched dynamics in solid 4 He. Consistent 
with earlier sections in this review, we assumed a Vogel-Fulcher-Tammann (VFT) form for 
the defect relaxation time t(T) = Toe A /( T ~ T °\ As in the previous sections, we obtain from 
our fits a negative 7b = —119 mK, in accordance with an avoided dynamic arrest of defect 
motion. 

A rough estimate of the electro-elastic coupling can be obtained from mass flow mea- 
surements in bulk solid 4 He, I2H1 where a pressure difference of A Pi ~ 0.1 bar across a 
centimeter-sized pressure cell was reported. The estimated local strain, with a bulk modulus 
B = 320 bar, is accordingly APi/B = 3 x 10~ 4 . This is consistent with the value we used for 
the fit in Fig.[T3] namely eo = 8.88 x 10~ 4 . In addition, the P{T) measurement by Yin clearly 
deviates from a purely Debye lattice behavior at around T = 0.4 K with a large positive in- 
tercept corresponding to the order of 100 ppm of TLS, s ee Fi g. [5] This number is roughly 
five times larger than the most disordered sample in L in's I ll2l] specific heat experiments on 
ultrapure 4 He with less than 1 ppb of 3 He impurities Ill2llll3 l. Thus the crystals grown by 
Yin are strongly disordered and harbor sufficiently many defects to support centers of local 
strain fields. 

Thus far, we assumed that both local and global stress are constant at low temper- 
atures. From Fig. [5] we can read off that the global pressure change between 300 mK 
and 40 mK is less than APj = 0.18 mbar. This is more than three orders of magnitude 
smaller than the local stress Oi = 8.88 x 10~ 4 x 320 bar ~ 250 mbar inferred from the 
dilatory strain eo used in the fit, as well as the static pressure difference APi =0.1 bar 
measured at two pressure gauges in mass flow experiments l28ll . Putting all the pieces to- 
gether, we find that the change of the dielectric function based on global density changes 
in the Clausius-Mossotti equation is negligible. This is because the corresponding den- 
sity change, Ap/p = APj /B < 10~ 6 , leads to a change in the dielectric function of only 
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8e « (e — \)Ap / p = Q.Q65APj /B < 10~ 7 , which is more than two orders of magnitude too 
small to account for the observed effect of order 10~ 5 . Whereas the model of local stress and 
electro-elastic coupling can account for the magnitud e an d temperature dependence of the 
dynamic dielectric function. More recently Yin et al. redesigned their dielectric func- 
tion experiment with a simplified capacitor geometry and found no measurable anomaly 
at low temperatures. Within our theory of disorder such a null result is consistent with a 
negligible amount of frozen-in stress in the solid. 



8 Conclusions and future directions 

In this review we provided a general overview of the role of defects in solid 4 He. We sug- 
gested that defect dynamics and freeze-out provides a rich ground to account for a significant 
fraction of the data, while at the same time they allow enough flexibility to accomodate sam- 
ple and history dependence and hysteretic behavior ubiquitously seen in experiments. The 
general response function approach presented covers a wide range of observed effects. We 
provide a brief synopsis of these below. 

8.1 Thermodynamics 

We started our discussion with the notion that any true phase transition, including super- 
solid, is accompanied by a thermodynamic signature. Therefore thermodynamic measure- 
ments like specific heat can reveal the signature one would naturally associate with a phase 
transition. However, we found that the excess specific heat and corresponding entropy are 
consistent with the contribution from noninteracting two level systems (TLSs). The esti- 
mated fraction of these TLS to the specific heat is at the level of tens to hundreds of parts- 
per-million. Consequently, the corresponding entropy contribution AS ~ 10~ 4 J/(K mol) is 
inconsistent by several orders of magnitude with reported values of the superfluid fraction 
(NCRIF) in the TO and shear modulus experiments. We also point out that some of the most 
disordered samples demonstrated "supersolid" fraction up to 20%, while there is little evi- 
dence for any excess entropy on the scale of the gas constant R times the fraction of defects 
in any measured sample. 



8.2 Single and double torsional oscillators 

We have shown that a phenomenological glass (viscoelastic) model for quenched defects 
accounts for the experimentally observed change in resonant frequency and the concomitant 
peak in dissipation. Our analysis of torsional oscillator (TO) experiments revealed that most 
are well described by a Cole-Cole distribution for relaxation times. In addition, we derived 
a simple relation for the ratio of change in dissipation and change in resonant frequency 
that can explain the large ratios observed in experiments. The values for the glass exponents 
in the distribution function of the backaction required to fit the experimental data point 
toward broad distributions of glassy relaxation times, which clearly invalidate any attempt 
to describe these experiments by a single overdamped mode, i.e., a single Debye relaxation 
process. 

We also applied these ideas to understand the double oscillator data. Our studies of the 
coupled oscillator showed that the observed shifts in resonant frequencies and dissipation 
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are in agreement with a glassy backaction contribution provided one includes anomalous 
damping in the dummy bob and an explicit frequency dependence in the backaction term. 
As a side comment, it came as a surprise that already the unloaded double TO (no 4 He) 
required a negative damping coefficient for the dummy bob to accurately describe resonant 
frequencies and dissipation at 300 mK. 

Finally, one should keep in mind that a significant difference between glassy and su- 
persolid dynamics is that a glassy contribution to the TO grows with frequency, while a 
superfluid component decreases with frequency. This could be another differentiating factor 
for separating very different relaxation mechanisms. 

8.3 Shear modulus 

We showed that the shear modulus anomaly of solid 4 He is strongly affected by the dynam- 
ics of defects. The freezing out of defects leads to stiffening of the solid concomitant with 
a peak in dissipation. By studying the glass susceptibility due to the backaction, we found 
that both the amplitude change and T-dependence of the shear modulus are well captured 
by this model. An important consequence of the dynamic response analysis was the descrip- 
tion of the dissipation or phase angle. We found that the peak height of the dissipation is 
independent of the applied frequency and linearly proportional to the Cole-Cole exponent 
a as well as the backaction strength go. As go depends on the concentration of the TLS, we 
predicted that increasing disorder will result in larger amplitude changes of the shear mod- 
ulus. Additionally, we extracted a universal scaling behavior proportional to cot(T) using 
the Cole-Cole plot. In this plot all curves of the shear modulus collapsed onto a single curve 
over a wide range of frequencies. Furthermore, we have shown that the glass contribution 
can be described by a viscoelastic model through the incorporation of anelastic elements in 
constitutive equations of stress and strain. 

8.4 Dielectric function 

We have shown that the arrested glass dynamics causes the low-temperature anomaly in 
strained solid 4 He through the acousto-optical (electro-elastic) coupling and proposed that 
the temperature behavior of the dielectric function is coupled to local strain fields near crys- 
tal defects. It records the glassy dynamics and freeze-out of the hypothesized TLS excita- 
tions, which also lead to a stiffening of the solid with decreasing temperature. This effect 
is not captured by the standard Clausius-Mossotti relation, which attributes dielectric func- 
tion changes to a change in mass density or polarizability of the nonpolar 4 He atom. An 
important consequence of the phenomenological glass susceptibility is the decrease of the 
dielectric function at high temperatures, accompanied by a broad dissipation peak that can 
be measured by the imaginary part of the dielectric function. We hypothesized that the coop- 
erative motion of atoms forming the TLS along dislocation segments is the relevant process 
contributing to the reported anomaly. In our model, both the change in e((B; T) and the dissi- 
pation peak are to leading order independent of the applied frequency. Since the coefficient 
go of the backaction depends on the concentration of defects, we predicted that the change 
in dielectric function will be larger in quench-cooled or shear-stressed samples, while it van- 
ishes in defect-free single crystals. Beyond the specific application to solid 4 He, which we 
invoked here, our formalism allowed for a direct link between elastic and dielectric proper- 
ties which could prove fruitful in many other systems. 
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In summary we encourage and welcome experiments that will allow a more precise 
structural characterization of 4 He. Also one needs to sharpen and provide a more detailed 
analysis of the structural aspects of solid 4 He in order to be able to investigate any ar- 
rested dynamics of the postulated glassy regions to separate it from a simple crossover 
phenomenon. Finally, we believe that more dynamic studies probing the frequency or time 
response to a stimulus are necessary to investigate the differences between small subsystems 
of glassy, supersolid or superglassy origin in solid 4 He. 
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